Spontaneous contractility— mediated cortical flow generates cell migration in 
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We present a generic model of cell motility generated by acto-myosin contraction of the cell cortex. 
We identify analytically dynamical instabilities of the cortex and show that they trigger spontaneous 
cortical flows which in turn can induce cell migration in 3-dimensional (3D) environments as well 
as bleb formation. This contractility-based mechanism, widely independent of actin treadmilling, 
appears as an alternative to the classical picture of lamellipodial motility on flat substrates. Theo- 
retical predictions are compared to experimental data of tumor cells migrating in 3D matrigel and 
suggest that this mechanism could be a general mode of cell migration in 3D environments. 

PACS numbers: 



Beyond its clear interest in the context of cell and 
molecular biology pQ , the study of cell motility and more 
generally the understanding of simple mechanisms of self- 
propelled motion is an important challenge for physics 
and biomimetic technology. Within this context, iden- 
tifying mechanisms of migration of micro-organisms and 
in particular living cells has motivated numerous works 
in the biology and physics communities [21 [3] . Sustained 
motion at low Reynolds number necessitates a constant 
energy input, and therefore an active system, that is a 
system driven out of equilibrium by an internal or an ex- 
ternal energy source. The cell cytoskeleton has been long 
identified as an example of such an active system. It is 
a network of semi-flexible filaments made up of protein 
subunits, interacting with other proteins such as motor 
proteins which use chemical energy of ATP hydrolysis to 
exert active stresses that deform the network pQ. 

Two principal out-of-equilibrium processes have been 
identified as responsible for cell motility: polymerization 
(treadmilling) of cytoskeleton filaments and contractility 
of the cytoskeleton, which results from the interaction be- 
tween actin filaments and motor proteins. Both polymer- 
ization induced motion [21 [3] and contractility induced 
spontaneous flows pQ have now been observed in-vitro, 
and studied theoretically [6] and numerically [7] . In 
all models the key ingredients of motion are an energy 
input to compensate dissipation and sufficient adhesion 
or friction with a substrate to acquire momentum. The 
usual picture of cell locomotion, following classical motil- 
ity experiments realized on 2-dimensional (2D) flat sub- 
strates, is then as follows: the cell lamellipodium builds 
strong adhesion points with the substrate and pushes its 
membrane forward by polymerizing actin. At the back, 
the cell body contracts and breaks the adhesion points 

mm- 

However, the geometry of the cell environment in vivo 
is significantly different from a flat substrate. The ef- 
fect of geometry - and in particular confinement - has 
been shown to play a crucial role in cell migration, en- 



abling the use of mechanisms widely different from lamel- 
lipodial motility [9l [10]. In this context, studying cell 
motility in 3D environments is a promising and widely 
unexplored field. In fact, recent observations [TT] reveal 
that MDA-MB-231 breast tumor cells migrate in 3D ma- 
trigel with a spherical shape (see Fig. [IJ according to 
a contraction-based mechanism in absence of actin pro- 
trusion or lamellipodium formation at the leading edge. 
Here, inspired by these experiments, we develop a generic 
model for motility based on the active contraction of the 
cell's actin cortex caused by the action of myosin. We 
analytically identify dynamical instabilities of the cortex 
and demonstrate that spontaneous cortical flows appear. 
We show that such cortical flows, observed in many con- 
texts of cell polarization and development [T21 [T3] , can 
also induce cell migration in 3D environments. Theo- 
retical predictions are compared to experimental data of 
|11) . and more generally suggest that this contractility- 
based mechanism could be a generic mode of migration in 
3D environments which strongly differs from the classical 
picture of lamellipodial motility. 

The model is as follows (see Fig. [I]). For the sake 
of simplicity we consider a cell of spherical shape and 
conserved volume with radius R, and use spherical co- 
ordinates (r, 9, 4>) (qualitatively similar results are ex- 
pected for other geometries). The cell cortex is a thin 
shell of acto-myosin gel which polymerizes at the mem- 
brane (which defines the outer boundary at r — R) with 
speed Up, and depolymerizes with a constant rate kd- 
Since the actin cortex can vary in thickness, it can be 
modeled effectively as a 2D shell of a compressible gel 
with density p(6,<f),t) — po(l + Sp(8,(f>,t)). The average 
density po is set by the polymerization speed v p accord- 
ing to po = aVp/kd where a is the actin monomer con- 
centration which for simplicity we assume is constant. 
Mass conservation is then given by Eq. ([!]) below. In a 
first approximation we write the pressure in the cortex 
P = aSp— f3W 2 Sp, where a -1 is the compressibility of the 
gel and j3 the correlation length of density fluctuations, 
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and denote v = (vg, v<f,) the components of the gel veloc- 
ity. Here, we consider only irrotational flows in the cortex 
and therefore write v = VV>- Hereafter the operator V 
implicitly acts on the angular variables 9, cf>. We neglect 
the effects of actin polarization and assume that the gel 
is in an isotropic phase. The activity of myosin motors 
then results in a diagonal active stress cr?g = —C,p5 a p 
in the gel which, following [14| , we assume is propor- 
tional to the local concentration of myosin in the cortex 
n(9,<p,t)i where £ is a phenomenological coupling. In 
the low Reynolds number regime, which is the relevant 
regime at the cell scale, the force balance given by Eq. ([2| 
is satisfied, where we denote by £ the friction with the ex- 
ternal medium. Here we assume for the sake of simplicity 
that viscous effects in the gel are much smaller than fric- 
tion and can be neglected as assumed in [T5]. The case 
of viscosity dominated dynamics, observed in |13j . can 
be studied along the same line and yields qualitatively 
similar results. In the cell bulk actin filaments are much 
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FIG. 1: Cell migration in 3D environments, a) Confocal sec- 
tion of an MDA-MB-231 cell migrating in 3D Matrigel. The 
cell expresses mCh-Lifeact, a fluorescent probe that labels F- 
actin White arrow indicates the direction of movement. 
Arrowhead points to the accumulation of F-actin at the rear 
of the cell. Blue arrows indicate the direction of observed cor- 
tical flow analyzed in Fig. [4] Scale bar, 10 pm. b) Cartoon 
of the model in the first unstable mode Yi,o, with the actin 
cortex in purple. Arrows indicate the velocity of the gel. 

more diluted than in the cortex and assumed to be in an 
isotropic phase, so that myosin, whose bulk concentra- 
tion is denoted by c(r, t), diffuses with diffusion constant 
D c (see Eq. (|3| where A denotes the 3D Laplacian). 
The myosin motors can attach to the cortex at rate fc on , 
where they are transported with the flow of the actin cor- 
tex itself and diffuse with diffusion constant Z? M , and can 
fall off from the cortex with rate k oS . Eq. Q describes 
the conservation of myosin at the cortex/cytoplasm in- 
terface, and Eq. ^ in the cortex. Finally, the non linear 
dynamical equations for the system are written: 



d t 5p + V ■ ((I + 5p)ViP) = -k d 5p 
V (-(> - aSp + PV 2 5 P ) = £V?A 
d t c = D c Ac 

- D c d r c\ r=R = k on c(R) - k oS n 



(1) 
(2) 
(3) 
(4) 



d t p + V • (pJVip) = k on c{R) - k oS p + Z^V V (5) 



The homogeneous static solution of Eq. (I|5| is triv- 
ially given by ip — 5p — 0, p — po and c = k^po/k 011 . 
In this case there is no actin flow in the cortex and 
the cell is at rest. To determine the linear stability of 
this homogeneous state we consider a perturbation of the 
(l,m) spherical harmonic for each variable. Namely we 
let c(r, (f>) = Co + ci :m (r)Yi ym (9,(j))e st where Co is the ho- 
mogeneous static solution given by Co = k° s p /k on , and 
define along the same lines the other (l,m) components 
of the perturbation as ipi,m, <5pz,m, Mz,m- Solving Eq. ([3| 
then gives c^ m (r) = c^ m I l+1 /z[Jif/)/y/r in terms of 
the Bessel function /„. We then obtain the following lin- 
earized homogeneous system: 



sdpi, m + k d opi, m — ipim = 

K 



(6) 



- CMi,m - CtSpijn - (3 - S pi !m - ^1pl,m = (7) 

D c ((/ + 1 - k on R/D c )I l+1/2 (u) - ttZi_ 1/2 («)) c hm (8) 
+ k oS R z l 2 pi, m = 

inff 1(1 + 1) v 1(1 + 1) 

{s + k otl + R2 ' D^)p ltm - fio R2 4>l,m 



(9) 



k on c lim I l+1/2 (u)/VR = 



where u = 



R. 



D Demanding a non zero solution 

yields an explicit equation for s which defines the dis- 
persion relation of the problem. This equation, which 
is non polynomial, can be studied numerically and fully 
characterizes the linear dynamics of the problem. A 
simplified (and analytically tractable) dispersion relation 
can be obtained in the limit where curvature effects can 
be neglected, which essentially amounts to considering 
the equivalent ID problem. This assumption is valid at 
time scales r 3> R 2 /D c , when the bulk concentration of 
myosin can be taken as constant. It can be shown that 
in this regime the dispersion relation is cubic in s and 
reads 



({s + D c k 2 ){s + k oS + D^k 2 



\s + D^k 2 



[!i{s + k d ) + k 2 (a + (3k 2 ) 

Cp k 2 (s + k d )( S + k on + D c k 2 ) = Q, (10) 

where fc 2 = 1(1 + 1) / R 2 . A numerical analysis shows 
that for the biologically relevant range of parameters (see 
figure [2]), this ID approximation is quantitatively very 
close to the exact dispersion relation, which means that 
in practice curvature effects can be neglected. This ID 



approximation ( 10 ) is useful to determine the threshold 



activity beyond which the instability appears. Indeed, a 
straightforward analysis shows that beyond the critical 
value defined by 



Cc 
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real positive solutions for s exist, which means that the 
homogeneous state is unstable and that cortical flows ap- 
pear. This threshold is controlled by both the diffusion 
properties of myosin, which clearly stabilize the concen- 
tration profile, and the friction which governs the ampli- 
tude of the cortical flow. Importantly, this threshold is 
independent of the physical properties of the gel, a and 
/?, which however enter in stabilizing terms determining 
the shape of the unstable pattern. Moreover, the anal- 
ysis of Eq. ( 10 ) reveals that s is maximized for a finite 
mode l max which will be dominant to linear order and will 
characterize the developing flow pattern at least at short 
times. Figure ^ shows that the most unstable mode 
^max grows with the activity. Note that a full non linear 
treatment, which goes beyond the scope of this paper, 
would be necessary to discuss the long time behavior of 
the problem. Interestingly, the specific limit of fast ex- 
change and fast depolymerization (k° n ,k oS ,kd large and 
D c = for simplicity) can be discussed at non linear or- 
der. In this case the dynamics takes the simplified form 



(12) 



which can be viewed as a modified Keller-Segel equation 
with drift — C/£V fJ- L 16J. In this limit it can be shown that 
the solution blows-up as soon as /i > — £-D M /C We expect 
that in the physical regime this singularity is damped 
leading to smoother asymptotic profiles. 
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FIG. 2: Dispersion relation s(l) in log scale: blue lines 
ID (where k = + 1)/R) and red crosses 3D. Pairs of 

curves in ascending order are for activities £ = —0.01, —0.05, 
— 0.1 kg/ims -2 . Inset : most unstable integer mode number 
^max (integer I when s(l) is maximal), against the strength 
of activity |£|. Other parameter values (order of magnitude 
estimates of biologically relevant values taken from [151 1171 - 
k on = 1 urns' 1 , k oii = 0.1s" 1 , D c = 10 /im 2 s _1 , 



D 



0.1s" 



1 /im 2 s 1 , kd 
P '= 1000kg/ims" 2 , £ = O.lkg^m" 3 ; 
/i = 10 4 /im~ 2 . 



1000 kgpm^s" 2 , 
, R = 10 /im and 



These results have several implications in the context 
of cell biology. Firstly, we show that such spontaneous 
cortical flows, which are observed in many contexts of cell 
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FIG. 3: Experimental data from [TT] for the intensity of 
mCherry-Lifeact labeled actin around the cortex of MDA- 
MB-231 tumor cells seeded in 3D Matrigel. Theoretical line 
in black shows the predicted unstable mode Y\fi{8). Red tri- 
angles show the average of 19 individual cells shown with grey 
crosses. The position 6 = is set at the maximum intensity 
along the perimeter and corresponds to the back of the mov- 
ing cells. 
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FIG. 4: Experimental data from [TT] for the velocity of actin 
around the cell cortex calculated from kymographs (position 
over time) of the data in figure [3] Red triangles show the 
average of 9 individual cells shown with grey crosses. The 
same reference 9 — as in figure [3] applies, with the black 
line showing the theoretical velocity of the predicted unstable 
mode Y h0 (8). 



polarization and development [12 [T3] , can induce motion 
when the cell is embedded in a 3D medium. Quantita- 
tively, the above analysis shows that for C > Cc the mode 
I = 1, for which /i, dp, ij) are proportional to Y\ a(ff), grows 
and creates a converging flow towards the pole 9 = of 
the cell where the cortex thickens, while the pole 9 = it 
depletes, as shown in figures ( |lppl ). In fact, the fric- 
tion law of Eq. |2| gives the tangential stress a nt on the 
sphere, which can be integrated to yield the force exerted 
by the medium on the cell in response to the spontaneous 
cortical flow. Assuming that the flow in the medium is 



4 



much slower than the cortical flow, this force is non zero 
and reads F = ^-h R 2 £vo, where ho is the cortex thick- 
ness and i>o is the amplitude of the cortical flow, showing 
that the instability of the mode Yi q induces cell motion. 
Note that the flow produced in the medium in response to 
cell motion induces a friction force —F, so that globally 
the cell acts on the fluid as a force dipole as expected. 

Secondly, we note that this cortical instability is closely 
related to bleb formation [T?J [2"3Tl25| . In the unstable 
regime, we find that the cortex is weakened, and could 
even lead to a hole, for instance at the pole 8 — it for 
the mode Yi.o(0)- Above a certain activity threshold, 
the pressure inside the cell will be enough to break a 
weakened cortex, forming a bleb. Larger modes I > I 
correspond to several simultaneous blebs. Our proposed 
mechanism of bleb formation is closely related to that of 
|14j . with the notable difference that their instabilities 
rely mainly on calcium dynamics. 

Finally we compare these theoretical predictions to 
experiments conducted on MDA-MB-231 human breast 
adenocarcinoma cells seeded in 3D matrigel. These cells 
were found to invade and migrate in 3D matrigel with a 
striking round morphology, and to display neither mem- 
brane extension nor blebbing at the cell front, in contrast 
with classical experiments on flat substrates. Interest- 
ingly, silencing the main actin polymerization nucleator 
found in lamcllipodial structures (Arp2/3) had no effect 
on the cell migration velocity, while it is known to be 
critical for migration on flat substrates. Moreover, inhi- 
bition of myosin II using blebbistatin almost completely 
abolished cell movement indicating that myosin activity 
is essential for movement. This shows that an alternative 
motility mechanism, widely independent of actin tread- 
milling and dependent on myosin contractility, must be 
involved to trigger migration. We suggest here that the 
cortical instability studied above, which relies mainly on 
cortical contractility could be the mechanism at work un- 
der these conditions. High resolution confocal imaging 
of live MDA-MB-231 cells in matrigel confirmed a highly 
non uniform distribution of cortical actin, which accu- 
mulates at the cell rear (see Fig{l]and pi]) 

in agreement 

with the model. More precisely, live imaging analysis of 
actin cortical structures indicated the quantity of cor- 
tical actin per perimeter length, quantified by p in the 
model. Figure ^ shows quantitatively the accumulation 
of actin at the back of the cell, and depletion at the front, 
in agreement with the predicted mode Yi Q which is un- 
stable for C > Cc- Analysis of kymographs of these films 
provided the flow of cortical actin. Figure Q shows that 
the velocity profile obtained is also compatible with the 
mode Yi,o, which as we have shown induces cell motion. 
In addition, blocking betal integrin function with anti- 
human betal integrin monoclonal antibody 4B4 resulted 
in a strong reduction in migration and in the presence 
of blebs. In our model this corresponds to decreasing 
the friction parameter £ and consequently the threshold 



value (c, which leads to a shift of the instability towards 
the higher modes Z > 1, corresponding to the formation 
of multiple simultaneous blebs. Together these experi- 
ments [IT] strongly suggest that in this example of cell 
migration in a 3D environment, motility is induced by 
the cortical instability described in the present paper. 

To conclude, we have presented a generic model of 
cell motility generated by acto-myosin contraction of the 
cell cortex. We identified analytically dynamical insta- 
bilities and showed that spontaneous cortical flow can 
induce cell migration in 3D environments. Theoretical 
predictions compare well to experimental data of tumor 
cells migrating in 3D matrigel, and overall suggest that 
this contractility-based mechanism of spontaneous corti- 
cal flows generation could be a general mode of cell mi- 
gration in 3D environments, which strongly differs from 
the classical picture of lamellipodial motility. 
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